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A direct derivation is given for the optimal mean fidelity of quantum state estimation of a d- 
dimensional unknown pure state with its N copies given as input, which was first obtained by 
M. Hayashi in terms of an infinite set of covariant positive operator valued measures (POVM's) 
and by Brufi and Macchiavello establishing a connection to optimal quantum cloning. An explicit 
condition for POVM measurement operators for optimal estimators is obtained, by which we con- 
struct optimal estimators with finite POVM using exact quadratures on a hypersphere. These finite 
optimal estimators are not generally universal, where universality means the fidelity is independent 
of input states. However, any optimal estimator with finite POVM for M(> N) copies is universal 
if it is used for TV copies as input. 
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I. INTRODUCTION 

iyy ' One of the essential differences between quantum theory and classical theories, from the information theoretical 
point of view, is that an unknown quantum state cannot be copied exactly, which was formulated as the no-cloning 
' theorem by Wootters and Zurek 0. 

Suppose we are given N identically prepared copies of an unknown state p on a d-dimensional space TL4 and try 
' to estimate the state p as precisely as possible by some measurement. Since we cannot increase the number of copies 
I by cloning the given unknown state, our performance surely depends on N, the number of copies given to us at the 
r — beginning. This is the problem of quantum state estimation, which has been studied since long ago 0, |j| . 

Two important points in formulating quantum state estimation are a prior distribution of the input states and a 
figure of merit to be optimized. In this paper we assume that the state p is pure and completely unknown in the sense 
that the state p is distributed over all pure states in a unitary invariant way. As a figure of merit we take the fidelity 
defined as tr [pp'\ — | ( (f> \ <// } | 2 , where p = \ (j> ) ( <j) | is the given input pure state and p' = \ (j)' ) ( <f>' | is the output pure 
f^S ■ state as a guess for p. 

""^ I The optimal mean fidelity in the case of qubit (d = 2) pure state estimation was found to be (N + 1)/(N + 2) by 
Massar and Popescu ■ They also pointed out that the optimal value of the mean fidelity is achieved by a joint 
measurement on the combined system of N copies but not realized by repeated separate measurements on each copy. 
Conceptually this unexpected result should be taken seriously since the input is a simple uncorrelated iV-fold tensor 
product p® N , though the improvement of joint measurement over separate measurement is relatively small (see also 
|5|). An algorithm for constructing an optimal and finite positive operator valued measure (POVM) has been given 
in 0. Bagan et al. also discussed optimal and finite POVM's in two-dimensional case with a different approach 0. 

Using the framework of covariant measurements |3|], Masahito Hayashi studied the estimation problem in more 
general settings: in general dimensions d and for a family of covariant error functions . He showed that the error is 
minimized by the unique infinite covariant set of POVM's in both Bayesian and minimax approaches, provided that 
the error function is a monotone increasing function of tr \pp'\ . As for the mean fidelity he found the optimal value 
to be (N + 1)/(N + d). Brufi and Macchiavello also obtained the optimal mean fidelity by establishing a connection 
between optimal state estimation and optimal quantum cloning [t| , the latter of which is another problem directly 
related to the no-cloning theorem. 

In optimal quantum cloning we are given N identically prepared copies of quantum state p on Jid and try to produce 
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a density matrix R p on H® M in an approximation of p® M as exactly as possible [lfj, [11|, [12J, [13| . There are two kinds 
of figures of merit for approximate cloning. In the many-particle test the full fidelity tr \_p® M R p ] is used for a figure 
of merit, whereas the one-particle reduced fidelity tr [pR p ] is employed in the single-particle test. 

The general formula for the optimal many-particle fidelity as a function of d, N, and M in the case of pure states 
was derived by Werner 0] It was also shown that the optimal fidelity is attained by the unique doner. This unique 
optimal doner was later shown to be also optimal with respect to the single-particle fidelity |13j . 

The connection established by Brufi and Macchiavello is the following (see also 11]). For given N copies of a 
pure state, first employ the optimal doner to produce infinite number of the best approximate copies, by which we 
can estimate the approximate copy as precisely as we want. On the other hand applying the optimal estimator to the 
input first, we obtain the best approximate estimation of the input by which we can produce infinitely many copies of 
the same quality. Thus they identified the optimal single particle fidelity in the large M limit with the optimal mean 
fidelity of quantum state estimation. 

For experimental implementation of POVM measurement, it is desirable that the number of outcomes of POVM 
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measurement is finite. However, finite optimal POVM's for state estimation have been constructed only in the two 
dimensional case (qubit) [EIS0E1- F° r general dimensions, the optimal fidelity was derived @,EJ, but finite optimal 
POVM's have not been discussed so far. In this paper we will show that one can construct finite optimal POVM's in 
general dimensions. We also show that the finite optimal POVM may be chosen to be universal, where universality 
means the fidelity is independent of input states. These are the main results in this paper. 

In Sec. II, we will first give a direct derivation of the optimal mean fidelity of quantum state estimation of a 
e?-dimensional unknown pure state with its N copies given as input. Our main concern is whether the optimal fidelity 
can be achieved by a finite set of POVM's. Therefore we do not assume the covariance of measurement, since the 
covariance implies an infinite set of POVM's, when input states are specified by a set of continuous parameters as in 
the case considered in this paper. We also avoid employing the optimal single-particle fidelity of approximate cloning, 
which is not straightforward to obtain. In Sec. Ill, we will study an explicit condition for POVM operators for optimal 
estimators. Using exact quadratures on a hypersphere, we establish the existence of a finite set of POVM's of optimal 
estimators. Covariant measurements are universal whereas measurements with finite POVM are generally not, where 
universality means the fidelity is independent of the input states. In this respect we will also clarify the conditions 
under which measurements with finite POVM are universal and therefore also optimal in the minimax problem. 

II. OPTIMAL MEAN FIDELITY 

Suppose we are given N identically prepared copies of a randomly selected pure state p = \ <f> ) ( <j> | on a d-dimensional 
complex Hilbert space Ttd and try to estimate the state p as precisely as possible by some POVM measurement {E a }^ =1 
on p® N . Since all inputs belong to the totally symmetric subspace of Hf N , the completeness relation of the POVM 

can be written as Y] x E a = Sn, where Sn is the projector onto this totally symmetric subspace. With the outcome 
of the measurement labeled with " a" , we infer that the state was a prespecified pure state p a = | <fi a ) ( (f) a \ . Our task 
is to maximize the following mean fidelity: 

A 

F(N,d)=Y / (tr[EaP® N ]trlp a p}), (1) 

a=l 

with respect to our strategy, the set of {E a , p a }a=i- I n the above equation < • • • > means an average over the input 
state p. 

We assume the input state is distributed over all pure states on TLd in a unitary invariant way. First let us fix an 
orthonormal basis {| i ), (i = 1, • • • , d)} in Tid and write a pure state | <f>) as \(f>) — Yl^—i Ci \ i }, where the coefficients 

Cj satisfy the normalization condition X)i=i c *i c i ~ 1- We assume that 2(i-dimensional real vector (3?Ci, ^sc^i—^...^ is 
uniformly distributed on (2d— l)-dimensional hypersphere. It is clear that the distribution defined above is independent 
of the reference basis. More precisely let {| i }, (i = 1, • • • , d)} be another orthonormal basis. Then for any function 
/, the following can be easily shown: 

As shown in the appendix, the average of a product of the same number of c's and c*'s is given by 

(ci,c* Ci 2 c* ■ ■ ■ Ci,c* ) = j— — ; — 4tt (sum of all contractions between i's and j's) . (3) 
ji j2 jii [d + l — ljl 

Using this formula Eq.Q and writing a density operator for a pure state as p — \4>)(4>\ — c i c j\ * )(j ; |; we 

obtain the following useful relation for the average of an iV-fold tensor product of identical pure density matrices: 

(P® N ) = ^, (4) 

where the sum of all permutation operators divided by AH is identified with Sn and djv is the dimension of the totally 
symmetric subspace, which is given by djv = tr [Sn] = N+d-iCd-i- 

It should be noted that the relation Eq.Q is a consequence of the unitary invariance of the distribution of p, 
which can be seen in the following way. For any unitary U on Hd we have (jj® N p® N U + ® N ') = (p® w ) , implying 
that the operator (p® N ) on the totally symmetric subspace of Ti.f N commutes with U® N for any U. Shur's lemma 
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then requires that (p® N ) be proportional to Sn , since JJ® N acts on the totally symmetric space irreducibly. The 
proportional coefficient turns out to be l/djv by a trace argument. Thus we obtain the formula of Eq.Q. 
Going back to the mean fidelity Eq.Q, we first rewrite it as 



A 



F(N, d) = £ (tr [E a p®( N +Vp a (N + 1)] ) , (5) 

a=l 

where the trace is taken over a total of N + 1 subsystems and the operator p a {N + 1) should be understood to act 
on the (N + l)th subsystem only; namely, for a single-particle operator f2 we use the following notation: fi(n) = 
2.®(n-i) ® O (g i®(N+i-n) _ Tjghrg the formula Eq.Q, we perform the integration over p to obtain 

1 A 

F(N,d) = - ytY[E a S N+lPa (N + l)}. (6) 

By tracing out the (N + l)th subsystem in the above equation, we finally obtain 



1 A \ ( N \ 

F(N,d) = S> ^ 1 + E^) 

v ' a— 1 L \ n— 1 / 



(7) 

where we used the following relation which holds for any single-particle operator 17: 



^S N (^ + J^(n)j 



tr N+ i [S N+1 n(N + 1)] = jf TT S\ ( 1 H- > ^ <>(//) | • is; 

Now it is easy to obtain an upper bound for the mean fidelity. Since p a {n) < 1, we have tr [E a p a {n)\ < tr [E a ]. 
Applying this inequality to Eq.(J7| and using the completeness of the POVM, we immediately find 

nN,d)<-^ = %±±. ( 9 ) 

dN+i N + d 

Equality in Eq.© holds if and only if tr [E a p a (n)] = tr [E a ] for n = 1, • • • , N, implying that E a is supported by the 
intersection of supports of p a (n), namely, E a is proportional to pf N . 
Let us write E a as 

E a = d N w a pf N , (10) 

where w a is a positive coefficient and the common factor d/v is introduced for later convenience. The completeness 
of the POVM implies dNj2a=i w aPa N — $n- Recalling the formula Eq.Q, we conclude that the necessary and 
sufficient condition for the POVM that achieves the upper bound Eq.@ is given by 



w aP r = { P ® N )- (ii) 



E 

a=l 

The right-hand side in this equation is the average of p® N defined as a continuous integration over a hypersphere, 
whereas the left-hand side is the sum of a finite number of sample density operators p® N with positive weights; 
namely, a continuous integration is replaced by a finite sum in Ea. (|llf> . which is a standard technique of numerical 
integrations (quadrature). Though a quadrature is in general an approximation, it may be exact for a certain class 
of functions. For example, the quadrature with a trapezoidal rule is exact for any linear functions. In this sense the 
condition (|ll|l means a quadrature on the hyper-sphere that is exact for p® N . Since quadratures with those properties 
exist as explicitly shown in the next section, we conclude that the optimal value of the mean fidelity is given by 

N + 1 

F optimai (N,d) = 1 ^. (12) 
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III. FINITE SET OF POVM AND UNIVERSALITY 



In this section we first show that we can construct a finite set of POVM's that achieves the upper bound of the mean 
fidelity Eq.@ or equivalently there exists a finite set {w a , p a }a=i which satisfies Ea. l(TT|) . We write \4> a ) — 2i=i c 1 1 * ) 
for p a = | <p a )( (f> a |. In terms of the expansion coefficient c", Eg. l|ll|l is equivalent to 

A 

E<« c "---« = MM---^i)' ( 13 ) 

Equation imposes no condition on products of different numbers of c's and c*'s. To make the subsequent argument 
simpler, however, we assume that they are zero like their exact average value. Then it suffices to show that there 
exists a quadrature with positive weights on the hypersphere S 2d ~ 1 that is exact for any polynomial of degree 2N . 
For a point x on S m , m = 2d, we write its polar coordinate parametrization as 

Xi = cos 6»i, 

Xi = sin 9 1 cos 02, 

X3 — sin 0i sin 02 cos 6*3, 

Xm-2 = sin 6»i sin 6*2 sin0 3 • • ■ cos0 m _ 2 , 
Xm-i = sin 0i sin 02 sin 6*3 • • ■ sin0 m _2 cos<ji>, 

Xm = sin 0i sin 2 sin0 3 • • ■ sin0 m _ 2 sin</>, (14) 

with the range of angle variables < 6i < 7r and < <j) < 2ix. The integration measure is given by the standard form: 

T (>TT flX I>27Z 

(20i sin" 1 " 2 0i / d6 2 sin" 1 " 3 2 • • • / tZ0 m _ 2 sm0 m _ 2 / dcj>. (15) 
Jo Jo Jo 

Let us consider the integral of x V \X2 ' ' ' Xm" with non-negative integers v K which add up to 2N . For each single 
integration we construct an n-point quadrature of the type J dxf(x) — Y^ K —i ^k/X^k) with positive weights uj k that 
is exact for the functions under consideration. We start with the 4> integration: 

dcj) cos Vm ~ 1 (j>sm l/m <j). (16) 



In the integrand we have the (v m -i + v m )th. power of e 1 ^ or e~ l< ^ at most. Therefore a simple trapezoidal rule, 
lu k = 2n/n and 4> K = 2ttk /n, gives exact results provided that v m -i + v m < 2N < n. We should remember that this 
integral vanishes unless both ^ m -i and v m are even. 
Next we consider the m _2 integral 

J m _ 2 = I d9 cos"— 2 0sin I/ — 1+,y ™ +1 0, (17) 
Jo 

where the subscript "m — 2" of the variable is omitted. Note that we can assume v m -i + v m is even since otherwise 
the whole integral is zero by the (f> integration alone, which is exact, whatever wrong results other integrations produce. 
When v m -i + v m is even, by setting cos0 = x we obtain 

1 

dxx^- 2 (l-x 2 )^- 1+ ^ )/2 . (18) 

-1 

Since the integrand is a polynomial of degree v m -2 + Vm-i + v-m, we can use, for example, the Gauss-Legendre 
quadrature formula of weights w% and points a;| (see [l4T | for example). In terms of variable the rule, u> K = wf / sin0 t; 
and K = cos -1 x|, gives the exact result provided that v m -2 + ^m-i + v m < 2N < 2n. Let us note that l m -2 = 
unless v m -2 is even. 

We must examine one more integral, the m -3 integral: 
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By the same reason as in the case of I m -2 we can assume that v m -2 + ^m-i + v m is even. Then the integration 
range can be enlarged to [0, 2ix\ so that an argument similar to that in the cj> integration applies. It turns out that 
the rule lo k — 7r jn and 9 K — 2tt(2k — l)/n is exact if v m -3 + ^m-2 + ^m-l + v m + 2 < 2(N + 1) < 2n. If we changed 
the integration variable by cos 6 — x, this rule would correspond to the Gauss-Tschebyscheff quadrature [j^j. The 
integral I rn -3 does not vanish only if v m -3 is even. 

It is clear that a similar argument also holds for remaining integrals; namely, for the integral I m -i the Gauss- 
Legendre (-Tschebyscheff) type of quadrature can be used when i is even (odd). The weights w a of the whole integral 
are positive since they are given by the product of weights lo's of the single integrations. Thus we can conclude that 
there exists a finite set {w a , Pa}a=i which satisfies Ea. ((TTJl . 

Now that we have shown there exists a finite set of POVM's that achieves the optimal mean fidelity of Ea. (|12[) . 
we study universality of finite optimal estimators. By universality we mean that the fidelity is independent of input 
states. The unaveraged fidelity of an optimal estimator for input p® N can be written as 



£ tr [E a p® N ] tr [p a p] = d N ^ Vatr [pf ("+ V<™>] • (20) 

a— 1 a— 1 

From this equation we find the fidelity is independent of p if and only if ^a=i w apa^ N+1 ^ = £/v+i/d/v+i, namely 

E^ +1) = (p» (W+I) ). (21) 

a=l 

This is a stronger condition than condition that is required for optimal estimators for N copies. Therefore optimal 
estimators are not generally universal, but any optimal estimators for M{> N) copies are universal if it is used to 
estimate N copies of an unknown state. 

A closely related question to this issue is the following. Suppose that for given N copies of an unknown pure state 
we first produce M(> N) copies by the optimal doner and then estimate the resulting state by the optimal estimator 
for M copies. What is the fidelity of this apparently detourlike two-step estimation procedure? Using the unique 
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optimal doner from N to M copies given by [lj , 

T(p® N ) = ^LS M (p® N <g> 1^-N)\ s (22) 
d M V > 

and an optimal estimator {E a , p a }a=i f° r M copies, we find that the unaveraged fidelity of the two-step estimation 
is optimal and universal: 

A A 

tr [E a T(p® N )] tr [p a p] = d N ^ w a tr [pf N p® N ] tr [p a p] - F optimal (iV, d). (23) 

a— 1 a— 1 



IV. CONCLUDING REMARKS 



In this paper we gave a direct derivation of the optimal mean fidelity of pure state estimation in a way we find 
simpler than the original ones In order to show the existence of the optimal estimator with a finite POVM, we 

avoided the assumption of covariance of the measurement and the use of connection to the optimal cloning fidelity. As 
a figure of merit we used the mean fidelity. It should be noted that the optimal fidelity is not changed if we take the 
infimum of fidelity for a figure of merit. As shown in the preceding section, in this minimax approach the condition 
for finite optimal estimators for N copies is Eq. I|21|) , which gives the universal fidelity, instead of Eq. . 

We showed how to construct a finite set of POVM's for optimal estimators by the use of an exact quadrature on a 
hypersphere. We expressed the integration on hypersphere as a multiple of single integrations, for each of which we 
constructed an exact quadrature rule. But this procedure does not generally give the minimal set of points on the 
hypersphere, or equivalently the minimal set of POVM's. In the case of qubit the minimal set of POVM's has been 
studied for several values of N [l5| . The weight w a obtained by our procedure depends on " a" . It may be desirable 
to have a constant weight since w a is equal to the probability of finding outcome " a" for the random input considered 
in this paper; w a =< tr [E a p® N ~\ >. Exact quadratures with a constant weight for polynomials of degree t on the 
hypersphere are called spherical i-designs. There is an existence result for all values of t in any dimension |l6j | , but 
explicit examples are in general not straightforward to construct. 
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APPENDIX A 

In this appendix we sketch a derivation of the formula Eq.@ based on a generating function for the readers 
convenience. In the text we considered the average of a function / of a normalized complex vector c — (c\, C2, • • • , Cd) 
and its conjugate c + : 

where 

/poo d 
dcdc+ = Y[d(Rci)d(%a)8(c + c-l). (A2) 

It is convenient to introduce a generating function G of A = (Ai, • • • , Xd) and its conjugate A + : 

G(A, A+) = J dcdc + e^ x+c+c+x \ (A3) 
so that the average ( c ii c ^ 1 c i 2 c ^ 2 ' ' ' Ci i c ji) ^ s calculated as 



I * * * \ 

\ c il c j 1 c i2 C j 2 ' ' ' °il c ji/ — 



9 9 9 9 -G(X,X + ) 



dXl dX h oxi dX 3l 



/G(0). (A4) 



A=0 



First we express the <5 function in the form of a Fourier transform as 5(c + c — 1) = j- J cL;e lu '^ c+c ^ . Then the 
integration over c and c + can be performed by a Gauss integral. The result is 

G(A,A+) = M! [ du —L e ^T£ke«<>«*), (A5) 

j A + A 

where e is a small positive constant, which should go to zero in the end. Expanding e and performing the u 

integration by a complex contour integral, we obtain 

00 

G(A, X+)=n d Y / ,, n , (-A + A)". (A6) 

^— ' n\(n + d — 1)! 

n— 

Now the formula of Eq.© can be shown as follows: 

/ , . .v (d - 1)! 8 d_ _d d_ + t 

^W'W /!(i + d-l)!flA? '"9A* 1 j 

— (sum of all contractions between Vs and j's) . (A7) 



(Z + d-1)! 

It is also easy to see that the average of any product of different number of c's and c*'s vanishes 
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